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Abstract 

The  self-magnetic  field  associated  with  the  current  in  a 
planar  diode  is  shown  to  reduce  electron  emission  below 
the  Child-Langmuir  current  density.  As  the  magnetic 
field  increases,  the  diode  current  is  limited  to  the  critical 
current.  Here,  a  1-D  analysis  is  carried  out  to  calculate 
the  suppressed  current  density  in  the  presence  of  a 
transverse  magnetic  field.  The  emitted  current  density  is 
found  to  decrease  modestly  until  the  magnetic  field 
associated  with  the  critical  current  is  reached,  at  which 
point  the  emission  shuts  off  abruptly.  The  1-D  analysis 
remains  valid  until  orbit  crossing  occurs  as  the  current 
approaches  the  critical  current.  The  minimum  diode 
length  required  to  reach  critical  current  is  also  derived. 


I.  INTRODUCTION 

High  power  vacuum  diodes  are  used  to  produce  intense 
electron  beams  for  many  applications.  When  self- 
magnetic  field  effects  are  negligible  and  the  cathode  is  a 
space-charge-limited  (SCL)  emitter,  Child-Langmuir  flow 
is  obtained[l].  Here,  the  effect  of  the  diode  self-magnetic 
field  on  the  locally  emitted  current  density  is  investigated. 
Results  presented  here  are  applicable  to  cylindrical  diodes 
with  an  anode-cathode  gap  D  that  is  much  less  than  the 
diode  cathode  radius  Rc.  In  this  case,  the  diode  can  be 
treated  as  planar  and,  when  the  magnetic  field  is 
negligible,  electrons  cross  the  diode  with  straight-line 
trajectories  as  illustrated  in  Fig.  la.  However,  when  the 
self-magnetic  field  associated  with  the  diode  current 
becomes  large  enough  to  significantly  deflect  the  electron 
path  as  it  crosses  the  diode  gap  as  illustrated  in  Fig.  lb, 
the  increased  electron  space  charge  density  in  the  gap 
above  the  emission  site  reduces  the  local  emission  below 
the  Child-Langmuir  current  density.  At  high  magnetic 
field,  the  diode  current  is  limited  to  the  critical  current 
Icrit,  which  is  obtained  when  the  electrons  reach  the  anode 
at  grazing  incidence  (or  equivalently,  the  “electron 
gyroradius”  equals  the  gap  size)[2,3].  Critical  current 
flow  is  illustrated  in  Fig.  lc. 

A  1-D  analysis  is  used  to  calculate  the  suppressed 
current  density  in  the  presence  of  a  transverse  magnetic 
field.  The  problem  is  similar  to  that  of  calculating  the 
Hull  current  in  a  crossed  field  gap[4].  In  treating  electron 
flow  in  crossed  field  gaps,  the  Llewellyn  approach  is 
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Figure  1.  Schematic  of  diode  electron  flow  for  (a)  Id0  « 

Icrit  ?  (b )  Id0  <:  Icrjt ,  and  (c)  Ido  lci;t . 


typically  used[4],  whereas,  here,  Poisson’s  equation  is 
solved  by  direct  integration.  Also,  in  the  case  considered 
here,  the  magnetic  field  is  produced  by  the  diode  current 
itself  and  increases  with  distance  along  the  gap  causing 
the  curvature  of  the  electron  trajectories  to  increase  with 
distance  as  well  (see  Fig.  lc). 

This  problem  is  of  interest  for  determining  the 
minimum  length  Lcrit  of  a  diode  designed  to  run  at  Icrit. 
This  situation  arises  when  an  electron  beam  pinch  is 
desired  to  concentrate  the  electron-beam  energy,  while 
minimizing  losses  to  ion  current  which  are  proportional  to 
the  diode  length  in  cylindrical  diodes[5]  (or  diode  radius 
in  pinched-beam  diodes[6]).  Although  this  analysis  does 
not  include  ion  current,  the  desired  minimum  diode  length 
with  ion  current  can  still  be  reasonably  estimated.  For  a 
nonrelativistic  planar  diode,  it  is  well  known  that  the  SCL 
electron  current  density  is  enhanced  by  a  factor  of  1.86 
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when  ion  current  flows[7],  while  the  scaling  with  diode 
voltage  and  gap  remains  the  same.  Thus,  it  is  assumed 
here  that  the  minimum  diode  length  for  bipolar  flow  L™”t 
is  to  first  order  simply  Lr^t  =  Lcrit/1.86. 


II.  SOLUTION  WITH  B0  *  0 


The  planar  diode  geometry  considered  here  is 
illustrated  in  Fig.  1.  A  voltage  V  is  applied  across  the 
electrodes,  and  a  uniform  magnetic  field  is  applied 
transverse  to  the  electron  flow  in  the  y  direction.  The 
current  ld0  is  emitted  over  an  axial  length  L  as  shown  in 
Fig.  la.  The  problem  can  be  cast  in  terms  of  the  electric 
potential  ®(x)  and  the  magnetic  vector  potential  Az(x)  =  - 
B0x,  where  B0  >  0,  0(0)  =  0,  and  ®(D)  =  V.  For  SCL 
flow  in  equilibrium,  Ex(0)  =  d®/dx(0)  =  0.  Because  y 
and  z  are  ignorable  coordinates,  their  canonical  momenta 
are  conserved.  An  expression  for  vx(x)  is  obtained  from 
conservation  of  energy.  Thus, 
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(i) 


vy(x)  =  0  ,  and  vz(x)  =  -©cex  , 

where  ©ce  =  eB0/mec,  e  and  me  are  the  electron  charge  and 
mass,  c  is  the  speed  of  light,  and  vx(0)  =  vy(0)  =  vz(0)  =  0. 
Note  that  for  B0  =  0,  vx(D)  =  v0  =  (2eV/me)1/2,  which  is 
the  speed  obtained  by  an  electron  in  crossing  the  gap. 
Because  vx2  >  0,  Eq.  (1)  shows  that  Id0  <  Icrit  (or  ©ce  < 
)  where 
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and  B0  =  2Id0/eRc  so  that  these  planar  results  can  be 
related  to  problems  in  cylindrical  geometry  with  Rc »  D. 

In  equilibrium  with  a  uniform  applied  magnetic  field 
(so  that  d/dz  =  0),  the  continuity  equation  shows  that  the 


Figure  2.  Plots  of  current  density  J0/Jcl  (solid)  and  x(1)/D 
(dash)  as  functions  of  ©ce/  ro™1  . 


current  density  Jx(x)  =  J0  is  constant  and 

ne(x)  =  — TT  ,  (3) 

evx(x) 

with  J0  >  0.  In  the  next  section  where  B0  is  allowed  to 
vary  in  z  along  the  diode  length,  it  is  assumed  that  the 
equations  for  vx  and  ne  in  Eqs.  (1)  and  (3)  still  apply  as 
long  as  d/dx  »  d/dz  ~  0.  Using  the  expression  for  ne 
given  in  Eq.  (3),  Poisson’s  equation  becomes 
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Flere,  however,  it  is  convenient  to  rewrite  the  equation  as 
d2S  _  8jteJ0 
dx~  meVs 

where  S(x)  =  vx2(x)  and  vx2(x)  is  given  in  Eq.  (1).  The 
first  integral  can  be  obtained  by  multiplying  both  sides  of 
Eq.  (5)  by  dS/dx  and  integrating.  This  yields 
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where  the  two  boundary  conditions,  S(0)  =  0  and  dS/dx(0) 
=  0,  have  been  used.  The  upper  sign  applies  in  the  case 
when  vx  is  increasing,  while  the  lower  sign  applies  when 
vx  is  decreasing.  This  occurs  when  Id0  — »  Icrit  where  the 
electron  orbit  has  been  significantly  deflected  by  the 
magnetic  field  so  that  energy  is  being  transferred  into 
axial  motion  as  the  electron  approaches  the  anode.  The 
second  integral  can  be  obtained  by  direct  integration, 
yielding  for  dS/dx  >  0 
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where  vx(x)  is  given  in  Eq.  (1)  and  S(0)  =  vx\0)  =  0  has 
been  used.  This  solution  applies  for  x  <  x(l1  where  x<n  = 
4ji2eJ0/meCBM’  is  defined  as  the  position  where  vx  stops 
increasing  (i.e.,  where  dS/dx  =  0).  The  value  of  x(1) 
depends  on  the  magnetic  field  strength  through  Jo(©ce) 
and  CBce  directly  where  ©ce  =  ©®  when  xU)  =  D.  It  can  be 
shown[8]  that  co^  =  n  ©“'/(tt2  +  2)1/2. 

When  xU)  >  D,  Eq.  (7)  applies  for  all  0  <  x  <  D.  When 

(i) 


<  D,  Eq.  (7)  applies  for  0  <  x  <  x1  and  Eq.  (6)  must 
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be  solved  again  (using  the  lower  sign)  for  x(1)  <  x  <  D. 
That  solution  is  then  matched  to  the  solution  in  Eq.  (7)  at 
x  =  x(1),  yielding 
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Assuming  that  x(1)>  D,  the  final  boundary  condition  ®(D) 
=  V  can  be  applied  to  Eq.  (7)  to  find  the  eigenvalue 


Jo(ff>ce)-  Using  J0,  Eq.  (7)  then  can  be  solved  numerically 
for  <t>(x).  When  x(1)  <  D,  Eq.  (8)  with  ®(D)  =  V  is  used  to 
find  Jo(ff>ce),  and  then  Eqs.  (7)  and  (8)  are  used  to  solve  for 
O(x). 

When  ff>ce  =  0,  it  can  be  shown[8]  that  the  classical 
Child-Langmuir  result  is  obtained  with 


J0(0)  =  JCL  =  - - 


MmJ  D- 

In  this  case,  nc(D)  =  n0  =  Jci/ev0.  When  coce  =  ca"“ ,  vx(D) 
=  0  and  the  electrons  reach  the  anode  at  grazing 
incidence.  In  this  case,  Eq.  (8)  can  easily  be  solved  for 
the  eigenvalue  J0(  co™' )  =  J0cnt,  yeilding 
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It  can  also  be  shown[8]  that  dJ0/da>ce  — >  -°o  as  roce  — >  a>"“  • 
Thus,  the  electron  emission  abruptly  shuts  off  at  o)ce  = 
co™1 .  The  solutions  of  Jo/Jcl  and  x<n/D  (for  x<1(  <  D)  as 
functions  of  coce/  ca"“  are  shown  in  Fig.  2.  Values  of  the 
various  parameters  are  also  shown  in  Table  1  for  the 
special  cases  a>ce  =  0,  ,  and  ro"“ . 


III.  MINIMUM  DIODE  LENGTH  FOR  Icrit 

If  laminar  flow  is  assumed  in  the  diode  gap,  then  the 
above  analysis  can  be  used  to  calculate  the  diode  current 
Id(z)  as  a  function  distance  along  the  cathode.  Id(z) 
increases  from  zero  at  z  =  0  to  the  total  diode  current  Id0  at 
z  =  L  as  shown  in  Fig.  1.  The  local  self-magnetic  field  is 
given  by  B0(z)  =  2Id(z)/cRc  and  determines  the  local 
current  density  from  Jo(roce)  as  given  in  Fig.  2.  Flere,  this 
local  current  density  is  denoted  by  J0(z)  and  dld(z)/dz  = 
2tiRcJo(z).  This  diode  configuration  is  illustrated  in  Fig. 
la  for  Id0  «  Icrit,  in  Fig-  lb  for  Id0  <  Icrit,  and  in  Fig.  lc  for 
Id0  =  Icrit-  The  minimum  diode  length  needed  to  obtain 
critical  current  then  is  determined  by 


Figure  3.  Plots  of  J0/Jcl  (solid),  Id/Icrit  (dash),  and  A/jiD 
(dot)  as  functions  of  z/L0crit  - 


Id(LCrit)  =  2jtRc|L“' J0(z)dz  =  Icrit  ,  (11) 

Jo 

where  Icr;t  is  defined  in  Eq.  (2).  A  zero  order  estimate  for 
Lcrit  can  be  obtained  by  assuming  that  J0(z)  is  uniform  in  z 
and  equal  to  JCL,  yielding  L0Crit=  9D/2p2,  where,  P  =  v0/c. 
Results  from  integrating  Eq.  (11)  are  plotted  in  Fig.  3  and 
show  that  Lent  ~  1.07L°crit-  This  provides  both  scaling  for 
and  a  reasonable  estimate  of  the  minimum  diode  length 
required  to  reach  critical  current. 

If  A(z)  is  the  axial  distance  that  an  electron  emitted  at  z 
is  deflected  by  the  magnetic  field  while  crossing  the  gap 
(see  Fig.  3), 

rT  rD  V  ( x) 

A(z)  s  Jo  vzdt  =  Jo  _7Tdx  -  (12) 

JO  JO  vx(x) 

where  T  is  the  gap  crossing  time  and  the  position  z  along 
the  cathode  is  associated  with  the  local  value  of  B0(z)  [or 
equivalently  a>ce(z)].  Again,  it  is  convenient  to  rewrite 
this  integral  in  terms  of  S  using 

dx  =  ±f327leJ°  Vs-4to2esl  dS  , 
l  me  ) 

from  Eq.  (6),  vx  =  S1 2  from  the  definition  of  S,  and  vz(x) 
=  -<Dcex  from  Eq.  (1)  where  Eqs.  (7)  and  (8)  are  needed  to 
write  x  in  terms  of  S.  When  xn)  >  D,  only  Eq.  (7)  is 
needed.  When  x,!|  <  D,  the  integral  must  be  separated 
into  two  parts,  using  Eq.  (7)  for  x  in  the  interval  0  <  x  < 
x11*  and  using  Eq.  (8)  in  the  interval  x(1)  <  x  <  D.  A(z)  is 
plotted  in  Fig.  3. 

The  analysis  presented  in  this  section  remains  valid 
until  orbit  crossing  occurs.  The  orbit  of  an  electron 
emitted  from  the  cathode  at  z  +  5  crosses  the  orbit  of  an 
electron  emitted  at  z  if  z  -  A(z)  >  (z  +  5)  -  A(z  +  5).  By 
Taylor  expanding  A(z  +  5)  for  5«  z,  it  can  be  shown  that 
the  condition  for  orbit  crossing  becomes  dA/dz  >1.  In 
order  to  evaluate  this  derivative,  expressions  for  dJ0/dz 
and  d(Bce/dz  are  needed.  First,  note  that  dJ0/dz  = 
(dJ0/dff>cc)da>ce/dz  where  the  solution  for  J0(a>ce)  is  shown 
in  Fig.  2  and  its  differentiation  with  respect  to  ©ce  is 
straight  forward.  Then  note  that  do)ce/dz  = 


Figure  4.  Plot  of  d(A/D)/d(z/L°crit)  as  a  function  of  z/L°crit- 


Table  1.  Special  values  for  various  parameters. 
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(2e/mce2Rc)dId/dz  =  4jieJ0(z)/mec2  where  Id(z)  and  J0(z) 
are  shown  in  Fig.  3.  A  plot  of  dA(z)/dz  is  shown  in  Fig.  4 
where  A(z)  has  been  normalized  to  D  and  z  has  been 
normalized  to  L°cr;t.  With  this  normalization,  the 
condition  for  orbit  crossing  becomes 

d(A(z)/D)>C=^T,  (13) 

d(z/L°cnt)  D  2p- 

Forthis  nonrelativistic  analysis  (i.e.,  P  «  1),  Eq.  (13)  and 
the  results  presented  in  Fig.  4  show  that  orbit  crossing 
only  occurs  very  close  to  Lcrit.  For  example,  orbit 
crossing  first  occurs  at  z  ~  1.05L°crit  for  p  =  0.5  (see 
dashed  line  in  Fig.  4)  and  rapidly  approaches  Lcrit  ~ 
1.07L°crit  (see  dotted  line  in  Fig.  4)  as  P  decreases.  Thus, 
this  analysis  is  valid  over  nearly  the  full  diode  length. 

As  discussed  in  Sec.  I,  the  minimum  diode  length 
required  to  reach  critical  current  when  ions  are  present 

decreases  to  L®  =  Lcrit/1.86  =  2.6D/p2.  This  assumes 
that  the  ion  current  is  much  less  than  the  electron  current 
and  that  D  «  R,  because  it  has  been  shown  recently  that 
the  ion  enhancement  factor  of  1.86  increases  with  aspect 
ratio  Rc/Ra  for  cylindrical  diodes  with  the  cathode  radius 
Rc  larger  than  the  anode  radius  Ra[9]. 


IV.  SUMMARY 

The  emitted  current  density  is  found  to  only  decrease 
modestly  with  magnetic  field  until  the  self-magnetic  field 
associated  with  the  critical  current  is  reached  where  the 
emission  shuts  off  abruptly.  The  1-D  analysis  remains 
valid  until  orbit  crossing  occurs  which  has  been  shown  to 
only  occur  very  close  to  the  full  diode  length  when  the 
current  approaches  the  critical  current.  Thus,  the  estimate 
obtained  for  the  diode  length  required  to  reach  critical 
current  is  accurate.  It  is  also  predicted  that,  when  ion 
current  flows,  the  minimum  diode  length  required  to 
obtain  critical  current  decreases  by  a  factor  of  1.86. 
Extending  these  results  to  the  relativistic  regime  and 


including  ions  in  the  analysis  will  require  a  numerical 
solution  of  Poisson’s  equation  and  is  the  subject  for  future 
work. 
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